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O^. Abstract 

^r\' We derive causal relativistic fluid dynamical equations from the relax- 

ation model of kinetic theory as in a procedure previously applied in 
the case of non- relativistic rarefied gases in Ref. [S] . By treating space 
and time on an equal footing and avoiding the iterative steps of the 

0\ ' conventional Chapman-Enskog — CE — method, we are able to derive 

causal equations in the first order of the expansion in terms of the mean 

CO ' flight time of the particles. This is in contrast to what is found using 

the CE approach. We illustrate the general results with the example 

C^^ ' of a gas of identical ultrarelativistic particles such as photons under 

^^ , the assumptions of homogeneity and isotropy. When we couple the 

fluid dynamical equations to Einstein's equation we find, in addition 
to the geometry-driven expanding solution of the FRW model, a sec- 
ond, matter-driven nonequilibrium solution to the equations. In only 

K> , the second solution, entropy is produced at a significant rate. 

1 Introduction 

Recent developments in the study of very hot objects, especially of the early 
universe, have underscored the need for an adequate description of relativis- 
tic fluid dynamics. Yet there are questions about the derivation of the basic 
equations from kinetic theory that are not completely settled. It may seem 
surprising that any doubt should arise in regard to fluid dynamical equa- 
tions whose provenance stretches back into the nineteenth century, but the 
application of fluid dynamics under extreme conditions may well require a 
reexamination of the foundations. And the answers to the questions raised 



may have interesting implications especially in the relativistic case. An ex- 
ample that has led to some discussion is the suggestion that bulk viscosity 
(the dissipative resistance to changes in volume) may give modify the ex- 
pansion rate in a self-gravitating gas [Ml EQl [311 [16]. We shall not delve 
into that controversial issue, but we will use the relativistic expansion of a 
self-gravitating gas to illustrate the content of the equations derived here. 
However, our main aim in this work is the derivation of causal equations for 
relativistic fluid dynamics. 

The failure of causality that was encountered in early derivations of the 
relativistic fluid equations from kinetic theory was foreshadowed in the non- 
relativistic version of those equations j 33[ [23] . Further indications that all 
was not well were revealed when measurements on the propagation of sound 
in rarefied gases failed to confirm the theoretical predictions based on the 
Navier-Stokes equations [37]. Those equations are commonly derived by 
the Chapman-Enskog procedure in the classical case [6]. The experiments 
helped to identify a weakness of the standard derivations of the fluid equa- 
tions [8] that make use of an iterative procedure whose ill effects are most 
pronounced when the mean flight times of particles are not sufficiently short. 
The approximation that we advocate here allows for the possibility of long 
mean flight times and the consequent deviations from equilibrium. 

Fluid dynamical equations may be derived as the the three lowest mo- 
ments over momentum space of the equations of kinetic theory, such as 
Boltzmann's equation. However, the stress tensor, which appears in those 
moment equations, is not prescribed and an additional equation for it is 
needed to complete the fluid description. A similar limitation occurs for the 
third moment of the kinetic equation — the heat equation — which contains 
no specification of the heat flux. An equation for that moment likewise must 
be prescribed before the system of continuum equations is complete. 

A standard way to prescribe the needed additional equations is the 
Chapman-Enskog approach in which an approximate solution to the kinetic 
equation is found by expanding the one-particle distribution function in the 
mean free path or mean flight time of the constituent particles. This ap- 
proximation is then used to obtain expressions for the pressure tensor and 
heat flux that may close the system of continuum equations. 

It has been the practice, following Chapman and Enskog, to simplify 
the derivations of expressions for the stress tensor at any order by intro- 
ducing results found in previous orders. In the nonrelativistic case this is 
done in such a way that, as Uhlenbeck [17] has described it, "the successive 
hydrodynamical equations are always of the first order in time derivatives of 
the macroscopic variables, but are successively of higher order in the space 
derivatives" (italics his). This leads to parabolic equations and that is at 
the origin of the failure of causality in both the classical and the relativistic 
cases [22] since space and time are not kept on the same footing in this 
approach. 



Grad [18] has observed that "the Chapman-Enskog expansion is asymp- 
totic rather than convergent." To us, Grad's suggestion imphes that the C-E 
iteration is not appropriate when the expansion parameter is not truly in- 
finitesimal. We have found support for this suggestion in the non-relativistic 
case where fluid equations derived from kinetic theory by omitting the iter- 
ative step of Chapman and Enskog give better agreement with experiment 
for the phase speed of ultrasound and of the structures of shocks than do 
the Navier-Stokes equations [U HH E] , when the mean free paths are long 
compared to macroscopic scales. 

Though not all difficulties have been resolved as yet [H], it appears worth 
treating the relativistic case without introducing the iterative approximation 
of Chapman and Enskog as we have done already for the case of radiative 
radiative fluid dynamics [9l[T0]. Here, we emphasize the formal structure of 
the derivations that lead to causal, covariant equations of fluid dynamics in 
the flrst order in the expansion in terms of the mean flight time of particles 
(to be denoted herein as e). In the relativistic case, many (if not most) 
derivations of the fluid dynamics equations from kinetic theory follow the 
procedure of Chapman and Enskog [47j . 

One result of applying the C-E procedure is that it leads in the first order 
of the expansions to non-causal equations that indicate (unphysical) insta- 
bility for the equilibrium state variables [331 ES EH ESI EH SO] • In previous 
work the difficulty has been circumvented by going to a second-order theory 
phenomenologically by Miiller [33] in the classical case and later extended 
to the relativistic case both phenomenologically and in kinetic theory by 
Israel and others in an ever developing literature beginning, for example, in 
[231 [2^ [36] . Also noteworthy is the version of Chapman-Enskoggery formu- 
lated by Van Kampen [49] who saw the development in terms of a singular 
perturbation theory based on the two-time approach. The meanings of fast 
and slow times that the scheme involves is not invariant under change of 
frame and the fluid equations found by Van Kampen are consequently not 
Lorentz covariant. 

In the derivation given below, we avoid the C-E iteration and so obtain 
relativistic fluid equations that are both covariant and causal in the flrst 
order of the development in mean flight time (or mean free path). For 
the purpose, we use the relativistic form of the relaxation model of kinetic 
theory, which is more general, if less explicit, than the relativistic form of 
Boltzmann's kinetic equation. Moreover, the derivations and their meaning 
are most easily appreciated in the simpler case of the relaxation model of 
kinetic theory. 



2 The Kinetic Model 

As our aim is to exhibit a procedure for deriving causal fluid equations, 
we shall consider only the relatively simple example of a gas of particles of 
a single kind having no internal degrees of freedom. Although the typical 
examples in which relativistic fluids are studied are astrophysical plasmas 
that are rich soups of particle species, we make this simplification in order 
to be able to carry out our procedure in a reasonably transparent manner. 
This approach may be appropriate when there is a dominant constituent in 
the mix of particles, and that is the view we adopt in the example discussed 
in the latter portion of this paper. 

We shall work (at first only) in Minkowski spacetime with coordinates rr'^ 
and adopt units in which the speed of light is unity. Each particle in our gas 
of identical particles has rest mass m (which may be zero). We describe the 
state of the gas with a scalar one-particle distribution function, f{x^,p^), 
where p^ is the four momentum of a particle. From the equation for / 
we wish to derive equations for the macroscopic fluid quantities including 
the local macroscopic fluid velocity, u^. Choosing a u^ is tantamount to 
specifying (what astronomers call) the local standard of rest, that is, a local 
rest frame. For ease of presentation, we postpone the prescription of u^ 
until we have developed the subject further. 

We describe the evolution of / by an equation of the form 

p^f,^ = a-Kf . (1) 

The quantity a is the rate at which particles are scattered into states with 
momentum p^ from some other momenta; in that sense, it is like a transi- 
tion probability. The second term on the right of ([1]) is the rate at which 
particles are being scattered out of the momentum state p^ with x being 
proportional to the inverse of the mean flight time of particles. The same 
equation has been used to study radiative transfer and, in the relativistic 
case, was discussed by L.H. Thomas [l5ll30l[l] who gave the transformation 
rules under frame changes for /, a and x. Boltzmann's equation, which is 
based on two-particle interactions, can be written in this form with a as a 
quadratic functional of / and x a linear functional of /. 

We presume the existence of a local equilibrium, f/Q\ , for which the right 
hand side of ([1])) vanishes. This property leads to the relation 

/(O) = ^ P) 

which is analogous to the Kirchhoff-Planck law of radiation theory. We 
adopt this expression for /(q\ with the understanding that it represents only 
a local equilibrium in which the macroscopic quantities it depends upon 
(pressure, temperature and so on) may themselves vary with position and 



time. Then ([T]) becomes the relaxation model 

P^/,M = x(/(0) - /) • (3) 

The nonrelativistic version of ([3]) was adopted for the study of material 
particles by Bhatnagar, Gross and Krook [5j and by Welander [53j. It is also 
known as the BGK model and, more recently, as the WKBG model. Krook 
had been a student of Eddington and he, at least, may well have had the 
radiative analogue in mind in thinking about this problem. In our work, we 
have benefitted from the formal equivalence of the various versions of the 
theory. In particular, we learn from Thomas' discussion of the relativistic 
transfer equation that k may be expressed in terms of its value in the local 
rest frame by the transformation rule 

X = u^p^k (4) 

where k is the inverse of the mean flight time evaluated in the local rest 
frame of the fluid. A hat placed on a symbol will denote evaluation of a 
quantity in the local rest frame throughout. 

The equation governing the relativistic relaxation model is then 

ep^/,M = n.p'^(/(o)-/) (5) 

where e = l/i<o • In our formulation of the equation, / is a scalar so that 
([5]) is manifestly covariant in keeping with the practice of much of modern 
relativistic radiative transfer theory [19]. In the present study, based on 
the relaxation model, we shall not allow for an explicit dependence of e on 
momentum and will take it to be a function of only the local macroscopic 
properties of the gas. 

In previous work with the relativistic relaxation model, Marie [32j wrote 
the relativistic relaxation model as in (I2D as if x were an invariant; he did 
not discuss the issue of its transformation properties. Though Anderson 
and Witting [21 [3] also eluded the issue of transformation properties, they 
did introduce the factor u^p'^ . Hence, if the 1/r in their treatment that 
apparently corresponds to that of Marie (and to our x) were replaced by 
the equivalent of x, their formulation would be covariant. 

3 Macroscopic Description 

In the Boltzmann version of the kinetic equation (with no quantum me- 
chanics), the right hand side describes the results of elastic, two-particle 
collisions. These conserve five quantities written in the relativistic case as 
■0"^, A = (0, 1, ...,4) = (/i,4), where -0^ = 1 and tjj^^ = p^. That is, multipli- 
cation of ([5]) by ifj followed by integration over momentum space gives zero 
on the right side and the resulting equations conserve the number, momenta 



and energies of the colliding particles. In the case of the relaxation model, 
we impose the same conservation conditions, by requiring that multiplica- 
tion by tp followed by integration over momentum space gives zero on the 
right hand side. Since we have assumed that e is independent of p^, this 
condition, known as the matching condition [SJ , may be written as 

J^^fdP = Ji;^f^o)dP. (6) 

Here dP = d^p/e is the invariant volume element in momentum space where 
d'^p is the three-dimensional volume element in momentum space [44 ^ [27] and 
e is the particle energy. (This result is discussed in terms of the geometrical 
basis of transport theory in [29].) 

The matching condition ([6]) leads to an osculating equilibrium distribu- 
tion function for which the fundamental invariants have the same values as 
in the actual state of the gas. We may then obtain the macroscopic equa- 
tions, by multiplying equation ([T|) by tp and integrating over momentum 
space. Before proceeding, we should recall that we are working just now 
in flat Minkowski spacetime. Hence we do not yet need to distinguish be- 
tween partial and covariant derivatives. To go to curved space or even to 
introduce curvilinear coordinates, we replace the partial derivatives by co- 
variant derivatives as appropriate and signify this by replacement of commas 
in the subscripts by semicolons. This replacement may be justified formally 
by including a term p^ df /dp^ on the left side of the kinetic equation and 
complementing it with the geodesic equation for p^ where the overdot de- 
notes differentiation with respect to proper time (or to an appropriate path 
parameter). The commas and semicolons then take care of themselves. 

When we multiply ([T|) by -0^ = 1, we find that 



p'^lf.dP = 0. (7) 

Since p^ is independent of x'^, this leads to 

m^ = (8) 

where 

ATM = [ pffdP . (9) 



Then, when we multiply by ■0'^ = p'^ and integrate, we find 

T^'',^^ = (10) 

where 

rp^^l^ = fp^p'^fdP . (11) 



It can be helpful to decompose these basic quantities. For that purpose, 
we employ the projection operator 

h^"" = gf"" - uf'u^ (12) 

where g^'^ is the metric tensor (for now, the Minkowski metric). The com- 
ponent of A^'^ along u^ may be identified as the number density of particles 
in the rest frame, 

N = u^N'' , (13) 

and its spacelike component (orthogonal to u^) is 

,p = hi'^N^ . (14) 

To separate T^'^ into its constituent parts we use the (Eckart [T7]) de- 
composition of p^ = (e, p) where p is the three-momentum and we write 
p = |p|. We let e = u^p'^ be the particle energy in the local rest frame 
of the fluid with a corresponding meaning for p. We further introduce the 
spacelike vector, Z'^, with the properties 

If^l^ = -1 and I'^Uf, = 0. (15) 

Then we may write 

pi" = eu^'+pl^' (16) 

and note that 

1% = "H^ = ^^ = -1. (17) 

pi pZ 

These considerations apply also for a gas of photons, in which case, p = e 
and m = 0. 

The stress tensor may then be written in component form as 



rpt^U ^ ^^^,^V ^ p^^^V ^ pU^I. ^ p^,U ^ 


(18) 


L*6 

E = 1 (?fdP F^' = f epl^'fdP P^" = j fm^fdP . 


(19) 


, of course. 




E = u^^u^T^^ F^' = h^"'uPnp P^"" = h^'PW"'Tp„ . 


(20) 



(Hats are omitted for cosmetic reasons.) Then ([8]) and (jlOp are written out 
as 

u^N^p + A^T? + J/;^ = (21) 

uP{Ev!'),p + EvJ'd + (F'^M^ + F^n'^),;, + P^^^^p = (22) 

where 'd = u^,p- A few more elementary steps are needed to complete this 
formulation of the basic equations, beginning with the specification of the 
macroscopic velocity. 



4 Frame Choices 



Equations (j2ip and (j22p provide relations among the macroscopic quantities 
N, E, u^, F^, J^ and P^'^ . We have yet to choose the frame in which to 
express these relations. That is, we need to specify u'^, which van Kampen 
|48j regards as a thermodynamic variable in its own right. And we would 
do well to choose a frame that simplifies the equations. Indeed, any five 
relations among the macroscopic variables specifies a frame, so many choices 
are possible, the two most widely considered being those proposed by Landau 
and Lifshitz [27J and by Eckart [17]. 

1. The Landau-Lifshitz Conditions 

Landau and Lifshitz specify a preferred frame by imposing these five 
relations among the macroscopic quantities: 

u^N^' = n^iVj) (23) 

and 

u.T'^'^ = u.Ti;^ (24) 

where the subscript (0) implies evaluation in local equilibrium; N'f'^-. 
and TJ^Q. are the appropriate moments of /(q) • 

In local equilibrium, to good approximation, the number current and 
the energy flux each vanish. Equations ()23p and (|24p then lead to 



(25) 



Thus, in going to the Landau-Lifshitz frame we transform away the 
energy flux and reduce (f22]) to 

u^iV_^ + iV?? + J^^ = (26) 

u^'{Eu'')^^ + Eu"^ + P^y^ = 0. (27) 

With ( )■ = u'^{ )_^ the second of these equations may be written as 
Eu^^'d + (Eu'^y + P'^y^ = 0. (28) 

It can then be broken down into components. 

On projecting ()28p in the direction of u^, we get the energy conserva- 
tion equation 

E + E^ + u^P^"^^ = . (29) 

When we project ([25]) with the operator /ip^ of ([ED, we obtain 

Eu^ + Pf""^ = , (30) 



N 


= 


^(0) 


E 


= 


^(0) 


rpti 


= 






ATj) and u^UyTi''' = 


= Uf.UuTl'^y Th( 




N = iV(o) 




E = E'(o) 




J^ = 0. 



which is the equation of motion. 

We shall use the Landau-Lifshitz frame choice here but it is also worth 
recording another much-favored choice, that of Eckart. 

2. The Eckart Matching Conditions 

In the frame chosen by Eckart, the matching conditions are given by 
N^ = Nfi. and u,,UyT^'^ = u,,UuT^^^. These conditions lead to 



(31) 



Hence, in the Eckart frame, the particle current is entirely convective. 
In this frame, the fluid dynamical equations are 

u^iV,^ + iVi? = (32) 

u^'iEu"),^ + Eu"^ + (F'^u" + F^u^"),^ + P^y^ = . (33) 

For either of these two choices of frame, as for any other one that is consistent 
with the physics, we thus have, to first order accuracy in e, causal, covariant 
fluid equations. However, these descriptions are incomplete until a means 
of computing P^'^ has been provided. A standard procedure for this is to 
try to relate P^'^ to the lower moments of /. Among the ways for finding 
that relation, the most direct is perhaps to derive an approximate solution 
to the kinetic equation for / and, from that, to compute the moments of 
/ up to P^"^ so as to close the system. In our procedure, we enlarge the 
usual notion of closure relations to allow a dependence on the derivatives of 
the lower moments as well. The appearance of derivatives in our first-order 
expression for the pressure tensor may be thought of as a process-dependent 
feature of our approximation. Our aim here is to illustrate this point and, 
with that intention, we next derive an expression for P^'^ in the case of a 
simple situation of astrophysical interest. 

5 Approximating / 

5.1 An Illustrative Example 

To illustrate our approach to the derivation of a closure relation without 
invoking the C-E iteration, we continue to consider the case of a gas consist- 
ing of one type of particle with no internal degrees of freedom. To simplify 
the presentation even further, we presume that the particles are ultrarela- 
tivistic so that their masses may be left out of account in a description that 
would be suitable for a photon-dominated medium. This approach leads 
to a closure approximation that is representative of more general situations 



without the necessity of performing the more arduous calculations that situ- 
ations with complicated particle mixes would call for. And even in the case 
of the photon gas, interactions are possible without the intervention of real 
massive particles, since very energetic photons may scatter off one another 
by creating virtual e"*"— e~ pairs. 

The commonality of this example with studies of radiative transfer in 
material media is also helpful in other ways, as we have already mentioned 
[45j . Just as for photons, we will not assume that the particle number is 
conserved (though we are not explicitly including any quantum mechanics). 
Hence we need consider simply 

r^;; = o, (34) 

to which our fluid dynamical equations boil down. Then, to complete the 
formulation, we complement equations ()29p and (|30p with an expression for 
P^^'^ in terms of E and E. 

To find a closure relation for P^^^ we introduce a series expansion of / 
in terms of e into ([5|) to find an approximate solution. Then we use that to 
evaluate and relate the relevant moments of /. For conventional radiative 
fluid dynamics, we have carried this procedure out in [9j. There we took 
u^^ as the velocity of the ambient medium and treated it as known. In the 
present example, there is no background medium and we are working with 
a gas of ultrarelativistic particles whose masses we neglect as for a photon 
gas. We now have in mind a u^ that is the appropriate velocity field of the 
fluid itself even in the case of a photon fluid. Although we then obtain a 
pressure tensor of the same form as in [9], the change in viewpoint makes 
it desirable to provide a sketch of the derivation since the closure relation 
obtained is central to the present worlo- 

A nice simplification of the ultrarelativistic problem is that a particle's 
energy is equal to the magnitude of its three-momentum. That is, e = p := 
IpI (with c = 1) where p^ = (e, p). Then, in thinking about the macroscopic 
aspects, it is useful to introduce the null vector n'^ such that n^e = p^ and, 
as in (fT6|) . to decompose it into 



n^ = u^' + 1^. (35) 

As before, 

/^//i = -1 and l^'u^ = . (36) 

Thus, as in (jlSp . the components of the stress tensor may be expressed in 
terms of the quantities defined in ([19]) . 



^Treating the photon gas as a fluid in its own right by working in the radiation's 
inherent reference frame is also a useful way to obtain approximate solutions in conven- 
tional radiative transfer theory. The relationship of that approach to the more familiar 
Eddington factor method will be described elsewhere. 
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5.2 The expansion of / 

To develop an expression for P^'^ by seeking an approximation for /, we 
return to the relaxation model ([5]). Since p^ = n^e and n^p^ = e, it may be 
written as 

en^/,/. = /(O) - / , (37) 

where, in the present instance, p^p^ = 0. We then seek approximate solu- 
tions to this equation in the form 

/ = /(0)+e/i + O(e')- (38) 

Though we shall not always indicate the O(e^) error terms, they are impor- 
tant and should be kept in mind. Then, from (|37|) . we see that 

/i = -^''/(O),/. • (39) 

With this approximation, we may evaluate the stress tensor but first we 
must specify /(g)- 

For the photon gas of the present example, the equilibrium distribution 

/(O) = [e^"'^- - 1]-' (40) 

is appropriate. This is the equilibrium solution for a gas of bosons with zero 
chemical potential. As usual, we use this expression in a local sense and 
allow that both /3 (= 1/T in suitable units) and u^ may depend weakly (to 
use this dangerous term) on x^. 

The relaxation equation is always dragging / toward the evolving equi- 
librium whose changes are connected to / itself through the matching con- 
dition, ([6]). The variation of /(q) with x^ is must be found and are not 
specified a priori; the changes take place through the dependence of /(q) on 
P and u^ each of which may itself depend upon x^. We take advantage of 
this structure of /(q) by using the chain rule to write 

5^/(0) = {T^^dT + e,^d,) /(o) (41) 

where dz means a partial derivative with respect to z. In higher orders, we 
would need to include derivatives with respect to n'^ and any nonconstancy 
of e would also make itself felt (just as in second order in the nonrelativistic 
case [S]). For particles with nonzero rest mass, we might need to include a 
chemical potential as well. 

A further simplification is that, as for the photon gas, /(g) depends on e 
and T only through the ratio e/T (Stefan's law) so that we may write 

5^/(0) = (r;.-r/%,^)9T/(o) • (42) 

And so we conclude that, to first order in e, 

/ = /(O) - en^ {T^. - TlPup,^) 9t/(o) + 0{^) . (43) 
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6 The Closure Approximation 

6.1 Zeroth order 

In leading order, with / = /(q) , we may use the last of (fT9]l to evaluate the 
pressure tensor, which becomes 

P^^- = jpHH'^f^.^dP. (44) 

Because the medium is locally isotropic in zeroth order, the angular integral 
in momentum space involves only products of the P. Such integrals are 
evaluated by expressing the results in terms of u^ and g^^'^ and making 
suitable choices for the coefficients. We find that 



j m"" 



dn = -^h'"' . (45) 



We then obtain (in zeroth order) the closure relation 

p^'^' = -Phi"' (46) 

where 

is identified as the pressure. On using the first of (fT9]) . we see that 

P=\e , (48) 

which is the usual equation of state for equilibrium radiation, where E is the 
energy density. Moreover, for the equilibrium ()40p we know that E = aT^, 
where a is the radiation constant, so that serves to tie in the temperature. 
We may also obtain the equations of motion for a perfect fluid at this order 
but we have no need of them and do not exhibit them here. 

6.1.1 First order 

Next we introduce (j43p into (J19p . Performing the angle integral is again 
straightforward and we find that 

I FrFrdn = — (hf^'hP^ + hf'Ph'"' + h'^'h"'') . (49) 

J 15 

So we obtain on using the various definitions that 

pi'^ = -Phf''' + /i^^/z (^ + -A+E'"' (50) 
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in which 

4/x 



5/.. = ^r^^^-^TXp,, (51) 

5 

is the viscous shear stress tensor where 

^f^upa ^ f^^,p^ua ^ j^f^af^up _ If^l^^ J^pa (52) 

O 

and /U = Pe is the viscosity. 

The middle term on the right side of ()50p does not appear in the Chapman- 
Enskog approach, which is applicable in only the case of very short mean 
free path. In that approximation, the system is very close to equilibrium 
and so (as we shall see in the example below) that extra term would be very 
small. But this term is central to the difference of our approximation from 
the usual first-order results; it is a process-dependent term that allows for 
deviations from equilibrium. 

To recover the first-order closure approximation found by using the 
Chapman-Enskog procedure, we need only neglect our extra term, which 
is of order higher than the first since /x is proportional to e. But that trun- 
cation leads to acausal equations as noted by Israel [22^ .23J who (as did 
others) preferred to go to second order to regain causality. 

Another way to avoid the causality problem is to use the moment method 
based on Grad's approximation |24^ I25j . but the large number of variables 
implicated in that procedure makes it hard to apply the results. Yet the 
key question is, how good are the results quantitatively in the various ap- 
proaches? For answer to this question, we have elsewhere turned to the 
nonrelativistic case for which empirical tests are available. These favor the 
method described here as reported in [H1I41]. for instance. 

6.2 The Fluid Equations 

To complete our statement of the basic equations, we note that the projec- 
tions of P^'^^p on u^ and hp^ are 

^^P% = ^-^{E - eQ) + E'^^^.up (53) 



and 



iV/* 



w 



V^^r. = — V(^ - ^Q) - ^^^{E - eQ),^ + hppE^^% (54) 



3 MPv -V. 3 

where 



Q = E + \e^ . (55) 



On substituting (j53|l and (j54|) into ()29j) and (|30|1 respectively, we find 
Q(l - U,) = -u^r.^^ (56) 
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and, since hfj_au" 



V(4S - eQ) = {E- eQ),^ - 3^5^^^,. (57) 



Since ()50p involves a derivative of E as well as E itself, it does not have 
the usual look of a closure relation — it is process-dependent. We may 
however eliminate E at the cost of some complexity. That is, we could 
replace E in (j50p by using (j29p . That would bring in u^P^'^u , which we 
could replace using ([53]) with ([55]) . That substitution in its turn brings back 
E but this time as a term O(e^) (since // = Pe). If repeated indefinitely, this 
substitution procedure brings in terms of all orders in powers of e. That is 
the reason that one may expect that ()50p can be a significant improvement 
over the standard closure formulae. However, we have preferred to leave 
well enough alone here and to retain the closed form with E rather than the 
infinite sum. This sequence of substitutions is different than the Chapman- 
Enskog iteration as there is here no truncation; in any case, we do not use 
it but merely mention it for clarification. 

7 Expansion of a Self-Gravitating Medium 

To provide a simple illustration of the content of the fluid equations we have 
just exhibited, we couple them to Einstein's equation and adopt the usual 
simplifications of elementary cosmology theory — homogeneity and isotropy. 
To avoid wandering too far into complicated (or unknown) physics, we keep 
to the model of ultrarelativistic, identical particles. Though this illustration 
is selected mainly for its lack of distracting details, it is so well studied that 
it makes for good comparisons with earlier work. We are in fact acting 
to some extent in the belief that sufficiently early on in our universe, all 
particles were ultrarelativistic. The model we have adopted for the basic 
fluid may not be too misleading in that case. 

7.1 Two Solutions 

We next study equations ([56]) and ([57|) together with the Einstein field 
equation 

C" = T'"' (58) 

where the units are such that the Einstein constant is unity. Under the 
conditions of homogeneity and isotropy we have a conveniently tractable 
example of flow described by five equations for the five unknowns E, H 
and u^. The configuration then resembles that of the simplest early models 
of cosmology. Comparison to the results from that topic will give a first 
indication of what modification the equations found here may introduce 
into the study of a self-gravitating medium in the case where the mean 
fiight time of the constituent particles need not be infinitesimal. 
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For an isotropic medium with no shear, we have H'"^ = and equations 
([56|) and (f57|) become 

Q(l -He) = (59) 

where H = t?/3 and 

ii^iAE - eQ) = {E- eQ),^ . (60) 

With R{t) as the scale factor, we have H = R/R. 

If we adopt the currently preferred cosmological model with zero curva- 
ture, the Einstein equation (j58]) reduces to [52] 

E = 3H^ (61) 

where E is the energy density. If there were any curvature, the extra term 
it would introduce would be significant for R close enough to zero. The 
rest-mass density has been left out of account since we consider only the 
case of a medium dominated by ultrarelativistic particles. 

Equations (j59p . (j60p and (j6ip . admit the following two solutions: 



1. Solution 1 

One solution to ([59|) is Q = 0. In that case (see ([55]) ). we have 

E + 4EH = , (62) 

which gives the familiar result that ER'^ is constant. If we assume 
that, as in equilibrium of a photon gas, E oc T^, we have T oc R~^, 
the well-known cooling law for equilibrium radiation in an expanding 
medium. With Q = 0, we conclude from ([60]) that u^ oc R~^ and we 
have recovered the features of the standard FRW solution. 

2. Solution 2 

The other possible solution of (|59p is 

He = 1. (63) 

This is really an approximate solution that is good up to errors of 
order e^. The nondimensional parameter, eH, is the ratio of the mean 
free flight time of the particles to the time scale of the macroscopic 
expansion and is a local Knudsen number in fluid dynamical terminol- 
ogy- 

Though it is not excluded that we may also have Q = in this sit- 
uation, we shall omit that case here for brevity. So we look only at 
the case where Q is not zero. The temperature is then not forced to 
behave in the way that it does in local equilibrium in the FRW model 
and it is undetermined as yet. 
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The properties of e now become significant. We have already assumed 
that e does not depend on p^^ exphcitly, though it may depend on 
local (in time) macroscopic properties of the medium such as temper- 
ature and pressure. Those variables, in their turns, will generally vary 
with R. Once all those dependences are specified, equation ([63|) is a 
differential equation for R. Its solution may be written formally as 

rR JO 

t = to+ e{T{R))^ (64) 

where Rq = R{to). If we specify the relation of e and T, ([M|l becomes 
an integro-differential equation for finding T{R). Then, we would need 
to solve dSSD and dSHD for R{t). 

We have not done that as yet as we do not know the nature of the 
interactions of the particles in the very early universe . Rather, if as 
is often done in using the relaxation model, we were to assume for the 
sake of discussion that e is constant, H would also be (nearly) constant 
in this solution. In that case, we would have 

R oc exp{Ht). (65) 

where we have taken to = ~oo and Rq = 0. Alternatively, we might 
have conditions in which solution 1 occurs right after f = so that 
there is a crossover to solution 2 at some t = to with appropriate 
matching conditions. In either case, the e- folding time of R in solution 
2 is the mean flight time of the particles. 

Though solution (j65]) has a certain appeal it does seem unphysical. 
In particular, according to ()6ip . constant H implies constant E. This 
awkward outcome originates with the assumption of constant e, which 
we have introduced only to get a rough idea of the possibilities of 
solution 2. But e = 1/k and k = Not where N is the rest frame 
number density of scatterers and a is the scattering cross-section per 
scatterer in the rest frame. In that case, we see that a more reasonable 
choice is to take a as constant. With N ex R~^, we then find that 

i^oct^, (66) 

where we have simplified this solution by choosing -R(O) = 0. 

These possibilities cannot be properly evaluated without a more explicit 
particle model that we shall not attempt to formulate here. The main point 
to be made is that the truncation introduced in the Chapman-Enskog proce- 
dure has filtered out one of the two solutions for the motion of our extended 
equations for the expansion of a relativistic, self-gravitating gas. And that 
second solution has a different physical origin than the usual ones. 
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Solution 1 may be considered geometrical in origin; after all, expansion 
can occur even in a cosmology without matter. Solution 2 however is matter- 
driven as we see from the close connection between the expansion rate and 
the mean flight time. In that sense, solution 2 is analogous to a shock wave 
whose thickness is typically within an order of magnitude of the mean free 
path of the particles in the fluid. But, in soution 2, this shock-like behavior 
takes place in time rather than in space so that the expansion time scale is 
analogously of the order of the mean flight time of the constituent particles. 
In the two situations, the macroscopic inhomogeneity (whether spatial or 
temporal) is determined by the microscopic behavior of the particles. A 
role of solution 2 in an expansion might be to provide a transition between 
two different behaviors of solution 1 when there is an interval in which a 
breakdown of equilibrium occurs. 

7.2 Entropy Production 

In solution 1, once the medium is in local thermodynamic equilibrium, the 
expansion will not destroy the equilibrium if the particles are ultrarela- 
tivistic. However, since the expansion rate and the mean flight time are 
comparable in solution 2, we cannot expect any form of thermodynamic 
equilibrium to be established in that case. Yet the relatively infrequent col- 
lisions that do occur will produce entropy. The only macroscopic effect in 
the present model that can account for this is the volume change which, in 
this case, is a disequilibrating effect through (what is known as) the volume 
viscosity (or bulk viscosity). In order to gain some idea whether solution 2 
may be of physical interest (if it may occur at all) we estimate its entropy 
production rate in order to form an impression of what its physical impact 
on the physical state of the gas may be. 

For the purposes of this brief look at dissipation occurring in solution 
2 for the expanding medium, we treat e as constant as we are after only 
a rough estimate here. In that case, as we saw already, the expansion is 
exponential with R oc exp{Ht) and H is roughly constant. If this phase of 
exponential expansion happens early on when H is large, equation (|6ip tells 
us that the temperature will then be high as well. 

To get some impression of the magnitude of such effects we use the 
formula for the entropy generation rate that has been computed in the study 
of radiative fluid dynamics, for instance in [S]. There, we find that 

S = ^^^ (67) 

with 

where S is the entropy. This is a single temperature result that resembles 
that found by Weinberg [52l |3] for a two-temperature medium; in each case, 
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^ vanishes in an FRW equilibrium state. That is, since i9 = 3H, the rate of 
entropy generation is 

S='-f(^H4j. (SO) 

which vanishes in solution 1. 

For this illustration, let us imagine that eH is small at the outset and 
that it grows slowly and passes through unity: the expansion evolves so as to 
make solution 2 possible. If a crossover to solution 2 does take place, there 
will occur a transition period in which the macroscopic evolution is driven 
by the interplay of the particle interactions and the large-scale expansion. 
It is internally consistent to take e, T and H as approximately constant 
before solution 2 gives way to solution 1. During the time interval when eH 
is of the order of unity, solution 2 would then lead to the (approximately 
constant) entropy generation rate 

A 4^EH , , 

5 = ^^. (70) 

In that case, we estimate the entropy at time t\ to be 

5 = (ti - to)5 , (71) 

where solution 2 applies from to to t\. To find the total comoving entropy we 
multiply by B^ [30] . Then, at the hypothetical instant at which we presume 
that solution 2 returns the baton to solution 1, we find that 

i?35 = /?3e3^(*i-*o)gti. (72) 

To get some idea of what kind of quantitative effect to expect we com- 
pute the entropy generated during the conventional estimate of the time 
interval sometimes adopted for allowing inflation to cure some cosmologi- 
cal ills. Thus we take the previously suggested [26] to = 10~^^s and ask 
what the entropy is at ti = 10~^^s. For this, we adopt for H and Tq the 
values 6 x lO'^^s"^ and lO^^K, respectively. We also have i?o = cto- The 
total comoving entropy generated in the interval to to t\ is then of the order 
V^^ J jK. This estimate is comparable to current estimates of the value of 
the entropy of our universe [26] and has also been found when the entropy 
generation is ascribed to matter/antimatter annihilation [3| or to bulk vis- 
cosity [30j in the standard solution. Of course, other times and values for 
the parameters could have been chosen; we have put values found in the 
literature into our estimates for comparison. If indeed solution 2 did arise 
in those very early times it may then have contributed something to the 
heating of the universe. 
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8 Summary 

Already in nonrelativistic kinetic theory, derivations of continuum equa- 
tions typically give rise to situations in which properties of the medium are 
propagated at infinite speed. Such undesirable results are produced by the 
unequal treatments of space and time in the Chapman-Enskog procedure 
|47| . When the same practice is followed in deriving continuum equations 
in the relativistic case, the C-E procedure again leads to acausal fluid dy- 
namical equations. Of course, the problem is more urgent in the relativistic 
case where the origin of the difficulty has been isolated by Israel \^22\ [25] and 
others. According to Maartens [30j, 

The problem arises from the first-order nature of the theory, i.e. 
it considers only first-order deviations from equilibrium. The 
neglected second-order terms are, in fact, a necessary require- 
ment to prevent non-causal and unstable behavior. These terms 
transform the equations governing dissipative quantities from the 
algebraic first-order type into differential evolution equations. A 
key feature of the second-order theory is that the equilibrium 
and dissipative variables are considered on the same footing, so 
that the theory is well suited to dealing with non-stationary pro- 
cesses, such as would occur in the early universe. 

The present work is based on the belief that the problem is caused by the 
iterative step taken in the C-E procedure, a step that was likely introduced 
to make it easier to derive continuum equations from Boltzmann's equation. 
The point made in the present work is that one should not make the iterative 
step of Chapman and Enskog in the derivation of the continuum equations. 
As the expansion of the one-particle distribution in mean flight time may 
well not be convergent [IE], the iteration that has become standard in such 
derivations is inappropriate. On omitting the iterative step of Chapman and 
Enskog, as we have done here, one obtains in first order an approximation 
for the one-particle distribution function that contains the derivatives that 
Maartens reports are lost in the C-E procedure. Formally, those terms are 
O(e^) in the present work, but it is acceptable to retain them since we are 
admitting to errors of that order in our first order approximation. 

In discussions of these issues one hears the sometimes held opinion that, 
in a first-order theory, all terms of higher order that appear should be sum- 
marily discarded. But that is not obligatory if the results contain the ad- 
mission that the error is formally 0{e^). This error signal does not make 
it desirable to introduce second order terms ad libidem, but neither should 
second-order terms that arise naturally through the formal development be 
discarded without a good reason. And that is the prescription we advocate. 

Of course, we have carried out the corrected procedure here for only the 
relaxation model of relativistic theory, but that is enough to make the case. 
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(The calculation for the Boltzmann equation is rather harder and we have so 
far performed it only for the nonrelativistic Boltzmann equation with results 
that will be reported elsewhere.) As was already seen in the nonrelativistic 
theory, eschewing the C-E iteration does not produce perfection but it does 
cure the worst problems coming from the C-E procedure as reported in 

[Hiiiiiiaiis]. 

By starting from the relatively simple relaxation model of kinetic the- 
ory, we have been able to bring out the essential difference between our ap- 
proach to the derivation of continuum equations and that of the traditional 
Chapman-Enskog method. Though the relaxation model is less specific than 
the Boltzmann equation as regards the nature of the collisions, it avoids the 
worrisome features of the relativistic collision process when the particles are 
not point masses. Moreover, it does not have the limitation of being confined 
to binary collisions. 

A possibly interesting sidelight that we have described here is that, when 
we introduce the Einstein equation together with the usual simplifications 
of elementary cosmological models, we find a second expanding solution in 
addition to the usual one of the FRW model for the case of ultrarelativistic 
particles such as photons. In this second solution, the time scale of the global 
expansion is of the order of the mean flight time of the constituent particles. 
This feature of the temporal behavior is somewhat reminiscent of classical 
shock waves whose spatial extent is comparable to the mean free paths of 
the particles of the medium. Analogously, solution 2 has an expansion rate 
comparable to the mean flight time in the particle interactions. 

In turning to cosmology for our illustrative example of an expanding 
medium we felt that, as Weinberg has put it [52], "the temptation to try to 
construct some sort of model of the very early universe is irresistible." In 
any case, whether there is any relevance to cosmology or not in our look at 
an expanding medium, the simplifications of elementary cosmology theory 
have at least provided us with a tractable illustrative example. 
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